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Abstract 

We propose a formula for the eigenvalue integral of the hermitian one matrix model with 
infinite well potential in terms of dressed twist fields of the su(2) level one WZW model. 
The expression holds for arbitrary matrix size n, and provides a suggestive interpretation 
for the monodromy properties of the matrix model correlators at finite n, as well as in 
the 1/n-expansion. 
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1 Introduction 



It has been known for some time that there is a close relation between the hermitian one matrix 
model and the conformal field theory of one free boson |2|)1 EH1 1221 1221 • O ne of the quantities that 
appears naturally in the matrix model context corresponds, in conformal field theory, to the operator 



h(a,b) = exp^Jj+(x)dx). (1.1) 

Here A is a complex number and we have expressed the free boson theory (at the self-dual radius) in 
terms of the su(2) level one WZW model with the standard basis of the weight one fields denoted by 
J + (z), J 3 (z) and J~(z). Since the operator product expansion of J + (z)J + (w) is regular, no normal 
ordering is required to ensure convergence of the integrals that occur upon expanding l|l.l|l . 

As will be explained in more detail in the next subsection, the matrix model analysis suggests that 
the operator I\(a, b) will exhibit interesting monodromy properties for the su(2) fields in the limit of 
large matrix size n. The aim of this paper is to elucidate these in terms of the above conformal field 
theory description. In particular, we shall propose that the operator l|l.lf) can also be expressed in 
terms of dressed twist fields (see section rOl below for more details) which will make these monodromy 
properties manifest. To leading order in l/n the relation between I\(a,b) and twist fields had been 
suggested before in |.'S-'SM 1 .'Sj : here we shall propose an exact relation for all n. 

1.1 Relation to matrix models and 2d gravity 

The matrix integral for the hermitian one matrix model is given by 

Z mm [t\ [n) = (const) /d$ e -£ teWr W = / dXt.-.dXn J] e-±Z n k=1 W(\ k ) _ (12) 

■' ' N i,3=l,i<3 

The first integral is over all hermitian nxn matrices, while the second integral amounts to expressing 
the first one in terms of the eigenvalues of <E>, and W(x) = J2 m >o t?nX m is the potential; more details 
can be found e.g. in the review The relation to the free boson conformal field theory is established 
by noting that the integrand of l|1.2fl can be written as a correlator of free boson vertex operators 
|2ElEniE21l221- Using the language of su{2)\ the expression is, 

(n\e- H J+(xi) • • • J + (x n )\0) = e~^^ w{Xk) J]^ - x 3 ) 2 , (1.3) 

i<j 

where H = — J c W (z) J 3 (z) ^ . This formula is most easily verified by employing the operator 

product expansion J 3 (z)J + (w) = (z—w)~ 1 J + (w) + 0(l) to commute e~ H to the right. Here we have 
assumed that W{x) is analytic on C: this allows us to use contour integrals, and guarantees that the 
only contributions come from the J + -insertions. Using charge conservation, the matrix integral can 
then be expressed in terms of the exponential Hl.lfl as 

A™ ( \ 

<n| e - ff / A (-co,co)|0) = ^-^Z mm [|] (n) . (1.4) 

The correspondence extends also to correlators. Let Z mm \t, /]^™- > stand for the the right hand side 
of 11.2(1 . with an additional factor / in the integrand. For the matrix model resolvent we then have 
the relation 



(n\e~ H J 3 (z)/ A (-oo,co)|0) = ^-^Z mm [i, tr^ - £W' \z)] (n) , (1.5) 



as can be verified by writing out the integrals explicitly and comparing the integrands. 

Provided the potential W(x) increases fast enough for x — > ±oo, the eigenvalues of the matrix $ 
will condense, in the large n limit, on one or more intervals on the real axis, the so-called cuts (for 



2 



details, consult e.g. The correlator Z mm \t, tr _* M 1 ^ has, again in the large n limit, square 

root branch cuts on these intervals. Equivalently, continuing J 3 (z) in (|1.5f) around an endpoint of 
a cut in the large n limit results in J 3 — * — J 3 . This motivates the idea to model the endpoints of 
matrix model cuts by insertions of twist fields cr(x) on the conformal field theory side 122] EH- Free 
boson twist fields have also been considered in the context of the integrable hierarchy approach to 
two-dimensional gravity in jlbl I2UI I2T] . 

We will investigate this effect in a simplified setting, in which we choose W(x) in l|1.2[l to be an 
infinite well potential, that is, W(x) — for x £ [a, b] and W(x) = +oo otherwise. This effectively 
restricts the eigenvalue integrations from the real axis to the interval [a, b] , so that the relation (|f .4fl 
takes the simpler form 

A" 
(27t)™ n! 

Matrix integration measures (or potentials) which force the eigenvalues to lie on a contour which 
has an endpoint on the complex plane, rather then to start and end at infinity, are referred to as 
ensembles with hard edges, see e.g. jSJ where such models are treated and more references can 
be found. From the conformal field theory point of view, the relation (|f .6fl is easier to analyse than 
(|f .4(1 because the potential term e~ H is absent. For example, for the one-point functions one finds 



(n\h(a,b)\0) = j^-^ Z^ . (1.6) 



(n|J 3 (,)/ A (a,6)|0) = ^ [tr _j_](n) 
(nillJ^/AMJlO) = Z^[det(z- M) ±2 ] (n) 



(1.7) 



It should be stressed that even if we choose the interval [a,b] to coincide with the location of a cut 
in a matrix model with analytic potential W(x), the large n expansion of the free energies 

ln((n\e- H J 3 (z)I x (a,b)\0)) and ln((n|e" H J 3 (z)I x (-oo, oo)|0» , (1.8) 

as well as their behaviour in the double scaling limit, will be different. The reason is that for subleading 
effects in n" 1 , the decay-behaviour of the eigenvalues just outside of the cut will be important, and 
that has precisely been cut off by the infinite well potential 1 . 

As we have explained above one may expect, given the monodromy properties of the matrix model 
correlators, that the operator I\(a, b) is proportional to the product of two twist fields to leading order 
in 1/n. For finite n, the monodromy properties are however quite different. Indeed, the correlator 
Z mm ||, det(z — M)]( n > is single valued for finite n (being a polynomial of degree n), but has jumps 
at the location of the cuts in the 1/n-expansion. This is an example of Stokes' phenomenon (for an 
exposition see for example [2j, |341 app. B]): the analytic continuation in z of an asymptotic expansion 
(the 1/n-expansion in this case) of a function f(n, z) can be different from the asymptotic expansion 
of the analytic continuation. 

The double scaling limit of the hermitian matrix model, with an appropriately tuned potential, 
describes a (p, 2)-minimal model coupled to Liouville gravity (see e.g. |11|1 or, equivalently, (p, 2)- 
minimal string theory (see 021 f° r a summary of recent developments). The above effect has been 
given an interesting target space interpretation in the context of minimal string theory |34| . There, 
the target space is identified with the moduli space of FZZT branes and it is shown that while in 
the semi-classical limit this moduli space becomes a branched covering of the complex plane, in the 
exact quantum description the moduli space is in fact much simpler, being just the complex plane 
itself. On the matrix model side, the analogue of the string partition function in the presence of a 
FZZT brane is the correlator of the exponentiated macroscopic loop operator, Z mm [t,dct(z — M)]^ 
pi l4()ll4"Tll3()| . As described above, this correlator is single valued for finite n but develops branch cuts 
in the large n limit. This is an example for how classical geometry emerges as an effective concept in 
string theory. 

As we will see, also the simplified matrix model with the infinite well potential that we consider 
in this paper exhibits this behaviour. The alternative formula for I\(a,b) in terms of dressed twist 
fields that we shall propose will then give a suggestive explanation of this phenomenon. 

x We thank B. Eynard for a discussion on this point. 
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1.2 Summary of results 



To explain more precisely our formula for I\(a,b) in terms of dressed twist fields, we first need to 
give a few definitions, starting with the relevant twist fields cr±\(z). Around an insertion of <J±\{z), 
the three sit(2)-currents have a Z2-monodromy 



(J + , J 3 , J~ 
If instead of J c one uses the basis 



) .— » (\- 2 J-,-J 3 ,\ 2 J+). 



K 



±(AJ+ + X- 1 ,;-) , JC^ = f (AJ + — A _1 J~) — J 3 ioiu = ±l 



(1.9) 



(1.10) 



the twist fields ct±a(^) have the standard monodromy K 3 i— > if 3 and if" i— > —K v . In particular, the 
field if 3 (z) is single valued, and the two twist fields a+\(z) and <j-\(z) are distinguished by their 
if 3 -charge (i.e. the eigenvalue of the K 3 zero mode), which is \ and — j, respectively. In the if-basis 
it is obvious that the monodromy around cr±\(z) amounts to an inner automorphism of su(2) of order 
two. 

Let us define an operator S\(a,b) which is very similar in spirit to (a special case of) the star- 
operators introduced in |38l l3~9"j . Explicitly, it is given in terms of twist fields and exponentiated 
J~-integrals as follows, 



5 A (a, b) = (b -a) * a+\(b) exp ( - 



1 

2^A 



J 



Ci 



(x)dx\ exp ( 



1 

2^A 



J-(x)dx]a-x(a)] . (1.11) 

/ J reg 



Here C\ is an integration contour from a to b passing above the interval [a, b], the contour C2 has the 
same endpoints, but passes below the interval, and [• • - ] re g refers to a prescription to regulate the 
first order pole in the operator product expansion of J~ and a±\ (see sections 15.21 and 15.41 below - ) . 
The contours are illustrated in the following figure, 

Im 



cr-x(a) 




(1.12) 



The dashed line represents the branch cut between the two twist fields. 

As we shall explain below in section I2"3l the product of the two twist fields a+\(b)<T-\(a) that 
appears in S\(a, b) can be expressed in terms of an exponentiated integral of the form 



a +x (b)a_ x (a) = (6-o)"* :ex p(^ (a J+(z) + A" 1 J-(z)) 



(1.13) 



Qualitatively speaking, the J - integrals in the formula (|l.lll) for S\(a,b) can be interpreted as re- 
moving the J~ part of this integral, leaving behind only the J + integrals that appear in I\(a,b). 
This observation motivates the following operator identity for the two rather different looking expo- 
nentiated integrals (f 1 . 1 1> and (|1 . 1 1|) - 



h(a,b) = S\(a,b) 



(1.14) 



This equality is the main result of our paper. We have no complete proof for it; the supporting 
evidence will be given in section |SJ 

As a consequence of (|1.14|l it is now possible to see that the operator I\(a, b) does indeed display 
the monodromy properties described in the previous section. Consider a correlator of the form 
(n\ (fields) S\ (a, b)\ 0), where (fields) stands for any product of the currents J c (z). The integration 
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ii) 



a 5 

(1.15) 

It turns out (see section I5~S|) that the part of the integral along the ellipse is suppressed by a factor 
r -2n £ or some r > o. This results in the approximation 

(fields)^ (a, 6) 1 0) = (n| (fields) 5^ unc (a, 6)|0) (l + 0{r' 2n )) , (1.16) 

where S^ runc (a, b) is defined as S\(a, b), but with the J~ -integrals taken only over the short horizontal 
contours shown in figure l|1.15l ii). 2 The difference between S\(a,b) and 5^ runc (a, b) in a correlator 
with an out-state of charge n is thus non-perturbative in 1/n. In particular, both correlators in 
(|1.16[) will have the same 1/n expansion (the correlators have to be normalised appropriately to 
allow a 1/n-expansion, see section I3T5|) . But since in S^ unc (a,b) the points a and b are no longer 
connected by J~ integrals, the monodromy of the currents J c (z) around the points a and b is just the 
Z2-monodromy (|1.9I) of the twist fields. On the other hand, the monodromy of I\(a,b) (and hence 
that of S\(a,b)) is not given by l|1.9jl . For example, in the presence of I\(a,b) the current J + (z) 
is single valued (since the operator product expansion J + (z)J + (w) is regular), while under (|1.9|) it 
changes to A~ 2 J~(z). In this sense, the 'correct' monodromy of the currents J c {z) in the presence of 
I\{a 1 b) is a non-perturbative effect and cannot be seen in the l/n-expansion. This is the same effect 
as observed in the previous section for the minimal string, albeit here in a different model, namely a 
matrix integral with hard edges. 

Further, in the large n limit itself, the ./"-integrals are suppressed altogether, and S\(a,b) can 
be replaced by a product of twist fields. Taking into account the need to regulate (|1.11|1 (see section 
15. 3D . we arrive in this way at the second important result of our paper, 

(n|(fields)S\(a, 6)]0> = n'^ 2^ e 3 C'(-i) (b-a)* (n|(nelds)o- +A (&) o-_ A (a)|0) (l + C^n" 1 )) , (1.17) 

where £ is the Riemann zeta- function. This shows that our formula for S\(a,b) also has the correct 
large n limit. 

In passing from I\(a,b) to S\(a,b) we have effectively decomposed the monodromy across the 
interval [a, b] into a product of three terms, one each associated with one of the lines in i|1.12fl (the 
explicit product can be found in l|5.2|) below). This is analogous to a method introduced in [TU| 
to analyse Riemann-Hilbert problems, which can also be applied to investigate the asymptotics of 
orthogonal polynomials (see @JlI|, and the lecture notes [H]). There, the orthogonal polynomials are 
encoded in a Riemann-Hilbert problem with an appropriate jump matrix across the real line, and 
their large-n behaviour can be found by manipulating the contour along which the jump condition is 
imposed in a way analogous to (|1.12l) . 

The paper is organised as follows. In section [21 we review how free boson vertex operators and 
twist fields can be expressed as exponentials of integrated currents. To compare the properties of 
I\{a, b) and twist fields, we calculate some correlators of sit(2)i-currents in the presence of two twist 
fields (section^ and in the presence of I\(a,b) (section^J). Finally, the definition of S\(a,b) is given 
in section [5] where also its properties are investigated. Section H3 contains our conclusions. We have 
also included two appendices where some of the more technical calculations are described. 

2 Representing fields as exponentiated integrals 

Let us begin by explaining how one can represent fields in terms of exponentiated integrals. While 
this may seem unfamiliar at first, there is at least one example where this construction is actually 
well known. This is the case of a free boson that we shall review first. 

2 As explained in section [5.3l in addition the regulator introduces an overall factor. 



contour (|1.12fl can be deformed as in figure i) below. 




truncate 
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2.1 The case of the tt(l) representation 

The free boson theory with field X(z,z) has an u(l) symmetry that is generated by a weight one 
current H (z) = i dX(z, z) with operator product expansion 

H(z)H(w) = - 1 . 2 +0(l). (2.1) 
2(z — w)^ 

In terms of modes, H(z) can be expanded as H(z) — 'Yl n H n z^ n ^ 1 . These modes then satisfy the 
commutation relations 

[H m , H n ] = - m S m ,-n ■ (2-2) 

The corresponding stress energy tensor is T(z) — :H(z)H(z): , where the colons denote normal 
ordering, i.e. 

: H{z)H{z) : = lim (h(w)H(z) - 1 ) . (2.3) 
«>— >2 V 2(w — z) z J 

The modes of the stress energy tensor T(z) = J2 n L n z~ n ~ 2 define a Virasoro algebra with c = 1; in 
terms of the modes H n we have 

L n = ^ : H m H n - rn : , (2.4) 

m 

where the colons denote here the usual normal ordering of modes. 

An (untwisted) highest weight representations of the u(l) theory is generated by a state that 
is annihilated by the modes H n with n > 0, and is an eigenvector of Ho with eigenvalue fi, 

H n \n) =/*<yn,o|M>, «>0- (2.5) 

The corresponding vertex operator will be denoted by , and can be described by the usual vertex 
operator construction 

V^z) = :e 2vX(z) :. (2.6) 

We would now like to express this operator in terms of the current H(z) of the conformal field theory. 
At least formally we can write iX(z) — J H(w)dw, and thus we should be able to write the vertex 
operator in terms of an exponentiated integral. However, the exponentiated integral will have a 
non- vanishing vacuum expectation value, and thus it will not just describe the field V^, but rather 
the pair of together with its conjugate V-^. Thus one is led to expect |2SI 

V„(b) V- tl (a) = (b-a)- 2 » 2 :cxp(2 A1 f H(z)dz) :, (2.7) 

where the prefactor is needed to produce the correct scaling behaviour, as will be discussed further 
below (see (I2.1()|) 'l. In fact, one can show that this identity holds in arbitrary correlation functions. 
(For a definition of conformal field theory in terms of correlation functions see for example |25|.) To 
this end one observes that the V± M satisfy indeed their defining relations (|2.5|l since one calculates, 
using Wick's Theorem, 



i=l 



n . b 

Y[H(ui)H(w) :exp(2ju / H(z)dz 

1 1 71 

( r ) ( TT H M ■■ ex p( 2 ^ / H(z) dz 

\w — b w — a/ \ V J a 



^— b -^^)\ll H ^ : ^{^ H{z)dz):) (2.8) 

i=l Ja 

1 n 1 



where we have used that 



r^ = -r-k-T— (2 - 9) 

a (w — z) z w — b w — a 
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Since the Ui are arbitrary, it follows that any correlation function of H{w) with the integrated 
exponential has only poles in (w — a) and (w — b) of order one; the fields at a and b are therefore 
highest weight states. It is also manifest from the above formula (by taking the contour integral 
around a or b) that their eigenvalues with respect to Hq are ±/i. 

The above analysis determines the exponential up to an overall function of b — a. This is fixed by 
considering the vacuum expectation value of Q2.7|l . which equals 

(V^b)V^(a)) = (b-a)- 2 » 2 . (2.10) 

This is of the form (b—a)~ 2h (as required by conformal symmetry) precisely for the choice of prefactor 
made in (|2.7|l . It is easy to see, by the same methods as above and using the definition of the stress 
energy tensor in terms of the u(l) field, that the fields that are defined by the right hand side of (|2.7|l 
also have the correct conformal weight. 

As a non-trivial consistency check, one can confirm that these highest weight fields then also give 
rise to the correct 4-point functions. To this end we consider 

(V^h) V-^ai) V„{b 2 ) V-„(a 2 )) 

2 / / r bi \ / r b2 \ \ ( 2 - n ) 

= (6i- ai )- 2 ^ {b 2 -a 2 )- 2v /:exp(2/x / H(z) dz): :exp(2z/ / H(w)dw) 
The correlator on the right hand side can now be easily evaluated and one obtains 



ETU = (61 - ai y^ (b 2 - « 2 r- 2 £ H|± 



11 

1=0 

(61 - ai)~ 2 ^ (62 - 02) " 



" J 

1=0 i=l Jai 




2f fu ,,^\-2v 2 (i b l ~ h 2) (Ol ~ 02) 



(ax - 62) [bi - a 2 ) 



2fAV 



(2.12) 



This then agrees with the known answer. 
2.2 Representations of su(2) 1 

At the self-dual radius the free boson theory is actually equivalent to an su(2) current theory with 
k = 1. The su{2) current symmetry is generated by three currents ^ and J 3 of conformal weight 
one with operator product expansion 3 

J + (z)J-(w) = + + 0(1) 

(z — w) z z — w 

J 3 (z)J ± (w) = ±L-W+0(l) (2.13) 

z — w 

j 3 (z)j 3 H = - k +0(1). 
2(z — wy 

The modes of J 3 and J ± then satisfy the commutation relations 

l^mi J-n\ = ^ J m + n + km # m ,_ n 

\J m ,J n ] = i J-m+ri (2-14) 
f/ 3 / 3 1 - -TO (5 



3 For fc = 1 this then agrees with the above u(l) theory by setting J 3 = H and = V±\. 
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At level k = 1, the su(2) theory has only two irreducible representations: the vacuum representation, 
and the representation with j = \- Here j denotes the spin of the su(2)k highest weight representation 
whose states of lowest conformal weight are labelled by \j, m) with m — —j, — j + 1, . . . , j — 1, j. For 
j = j there are only two highest weight states, namely \\,\) and ||, — h) that are conjugate to one 
another. It is easy to see that the corresponding vertex operators V(\j,m),z) can then be obtained 
by the previous construction provided we take /i = ±^ 

V(\H),b) F(|!,-!), a ) = (&-a)-5 :exp( fj 3 (z)dz) : . (2.15) 

For larger values of k > 2, however, the above construction does not account for all non-trivial 
representations. In fact, in addition to the vacuum representation (that corresponds to [i = 0) only 
the highest weight state ||, |) together with its conjugate state || , — |) can be described in the above 
manner (with \i = i). 

2.3 Twisted representations of the sw(2)i theory 

Up to now we have only discussed the untwisted highest weight representations of su(2). The corre- 
sponding vertex operators have the property that the currents J 3 and J are single-valued around 
the insertion point of the vertex operator. As is well known, the su(2) theory (like any affine the- 
ory) has also twisted representations for which the currents have non-trivial monodromies around the 
insertion points of the vertex operators. For su(2) all of these twisted representations are actually 
equivalent (as affine representations) to untwisted representation, since all automorphisms of SU(2) 
are inner. However since the relevant identification modifies the energy momentum tensor, twisted 
representations describe often different physical systems (for an introduction to these matters see for 
example [26| ) . 

In the following we shall mainly be interested in Z2-twisted representations of su(2)i. One class 
of Z2-twisted representations have the property that the monodromy of the currents is described by 
the (inner) automorphism 

J 3 ^J 3 , J±h^._J±. (2.16) 

The corresponding representation then has modes J 3 that are integer valued (n G Z), while the 
modes are half- integer valued (r£Z+|); these modes then still satisfy the same commutation 
relations (|2.14|l as above. Since these twisted representations are in one-to-one correspondence with 
the usual untwisted representations, there are two inequivalent irreducible Z2-twisted representations 
for su(2)i: they are generated from highest weight states a± with Jq(t± — ±4<t±. From what was 
explained above, it is then clear that these vertex operators can also be described by exponentiated 
integrals; indeed we have simply 

a + (b) a_(a) = (b-a)-i : cxpQ J J 3 (z)dz): . (2.17) 

In the following another class of Z2-twisted representations will play an important role: these are 
the Z 2 -twisted representations for which the monodromy is described by 4 

J 3 ^-J 3 , ji^jT, (2.18) 

Obviously, this only differs by a field redefinition from (|2.16|) : if we define (cf. i|1.10[l with A = 1) 

X 3 = i(j+ + J-), K ± = ±i(j + - J~) — J 3 , (2.19) 

then the fields and K 3 satisfy the same operator product expansion as and J 3 (and thus their 
modes have the same commutation relations as l|2.14|l '). Furthermore, the monodromy l|2.18|l has the 

sec Ot with A = 1. The case of general A differs from A = 1 by the rescaling t—f A^ 1 J^. For the following it 
is therefore sufficient to consider A = 1. 
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same form as ()2.16[) . In particular, the two highest weight states a± that are now characterised by 
the condition that Kq<j± = ±j<j± are given by 

cr+(6) cr_(a) = (b-a)-i : exp( j j (J + (z) + J~{zj) dz\ : . (2.20) 
This is the formula that will motivate our ansatz for Sxidjb) (see section 5.1). 

3 Properties of the sw(2)-twist fields 

As explained in the introduction, we propose that the operator I\ (a, b) given in can be expressed 
in terms of twist fields as in To support this claim, we shall compare, in section^ a number of 

properties of I\ (a, b) and S\(a, b). As a preparation, we now want to study the correlation functions 
of the su(2) twist fields. 

3.1 Monodromy and zero-point function 

We will start with the zero-point function (n\<j + (b)(j_(a)\0) . Here (n\ denotes an out-state which is 
highest weight with respect to J 3 and has J 3 -charge n, 

(n\ = n5 m , (n\ for m < 0. (3.1) 

For n > one can write (n\ explicitly in terms of modes (for n < one has to use J + ) 

H=(0|JrJ 3 ----J 2 -_ 1 . (3.2) 

Let us normalise (0|0) = 1. One verifies that with the above definition (n\n) = 1, as well as 

(n\ J~ =0 for m < 2n and ( n \^m — f° r m — ~ 2n , (3-3) 

where we have used the null-vector relations of sm(2)i. The highest weight property (n\ Jf n = for 
m < and {n\jQ = n{n\ are then immediate consequences of the commutation relations (12.14(1 . In 
order to evaluate (n|<7+(&)cr_(a)|0) one can write the out-state in the form l|3.2fl and express the J~ 
modes in terms of the i-C-basis (|2.19|l . For example, abbreviating \aa) = a+(b)a-(a)\0), 

(l\aa) - (0|Jr|<r<7) = (0\(Kf - \{K+ - R-))\aa) . (3.4) 

As is clear from the discussion in section 12731 the field K 3 has only a simple pole with a±, 

K 3 {z)a ± {a)=±\j-^- )+ 0{l), (3.5) 

and thus it is easy to evaluate the term involving giving (0\K 3 \<j<t) = \{b- a)(0\aa). For K± 
one can write (c/. also ^3) 

(0|Kf| ff «7) = j c z(Q\K±{z)\*a) — = ^ yj{z -a)(z- b) (0\K±(z)\vv) — = . (3.6) 

To see the second equality, expand the square root in z" 1 and use the highest weight property of (0|. 
In the third step the contour Coo is deformed around the insertions tr+ (6) and <r_ (a) and the highest 
weight property of the latter is used. 

Altogether we thus obtain (1|<t<j) = j(b — a){0\aa). For higher values of n the calculation is 
similar, but more tedious. Luckily, the exact answer is known from factorising the correlator of four 
twist fields OES], [T71 eqn. (4.17)], 

(nK(6V_(a)|0)=4-« 2 (6-a)" 2 -S, (3.7) 

_ i 

where we normalised the twist fields such that (0|crer) = 1 • (b — a) 8 . 
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3.2 Correlators involving sw(2)-currents 

The next correlator we consider is the one-point function of J 3 {z) in the presence of two twist fields. 
It is easily determined by considering the function 

f(z) = ^(z-a)(z-b) (n| J 3 (z)<r + (b)a-(a) |0> (3.8) 

and noting that f(z) is single valued on the complex plane and does not have any poles. Since 
(n\jQ — n(n\ we furthermore have lirn^oo f(z) = n (n\ a + (b)<j-(a) |0) so that f(z) is in fact a 
constant. In this way we find 

(n\J 3 (z)a + (b)a.(a)\0) = " (n\ <r+(b)a_(a) jO) . (3.9) 

y?(z ~ a)(z - b) 

Correlators with several J 3 insertions can be determined in a similar fashion. 

Finally we will need correlators with several J insertions in the presence of two twist fields. 
These can be determined from the Knizhnik-Zamolodchikov equation, that is, by solving the first 
order differential equations resulting from the null-vector 

{L- 1 -2vJ 3 _ 1 )\J v ) =0, where v S {±1} . (3.10) 

Knizhnik-Zamolodchikov equations in the presence of twist fields have been studied in [5]. To find 
the differential equations first note that as a consequence of this null-vector we have the identity 



J (u—a)(u—b) T , . . TI . , . du vt 2z - a - b r — — d \ TI/ , , ,„ , , . 

^ '- J 3 (u)r(z) ^- = ^[ — r. — — + y/(z-a)(z-b) j- \j u {z), 3.11 

u-z 2-ni 2\^/{z-a)(z-b) 9z I 



where C z is a contour winding closely around the point z. Here we have used the operator product 
expansion 

J 3 (u) J v (z) - — ^— J v (z) + U 3 J v )(z) + 0(u - z) . (3.12) 

u — z 

Around a point w ^ z and around infinity we find analogously 



\/(u-a)(u-b) 3 du ^{w-a){w~b) 

J [u)J (w) - — = v J (w) 

u-z y ' K ' 2ni w- z y ' 

y/(u-a)(u-b) 3 du 

\ n \ J ( u )tt- = n(n\. 



(3.13) 



u — z 2ni 

Consider now the integral 



y/(u-a)(u-b) 3 du 

{n\J°(u)J x (zi) ■ ■ ■ J m (z m )\aa) — . (3.14) 

c ^ u — Zfe Im 

This contour integral can be calculated in two ways: on the one hand, we can directly use H3.11|l and 
thus evaluate the contour integral in terms of the right hand side of 13. lH . On the other hand, we 
can deform the contour around Zk to encircle all other insertion points Zi, i ^ k as well as infinity 
and a and b. As in the previous calculation, there is no contribution from the twist-field insertions at 
a and b. The individual contributions from the points z% and infinity can be evaluated using (|3.13(l . 
and one thus arrives at the system of partial differential equations 

D k {n\r i {z l )---J Vm {z m )\(j(j)=Q for k = l,...,m, (3.15) 

with 



v k ( 2z k -a-b r- 3 ^ ^{zi-a)(zi-b) 

Dk = T\ Ti \i IT + V (zk-a)(z k -b) — - n + }^ vt . (3.16) 

2 V V \ z k-a ){Zk-b) az k z t - z k 
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To solve the equations ((3.15(1 it is convenient to pass from the complex plane with twist-field insertions 
at a and b, which generate a branch cut on the interval [a, b], to a double cover (c/. The 
letters z, w will always refer to points on the complex plane which forms the base of the double cover 
and £, £ to points on the double cover. Since we have a single cut, the double cover has again the 
topology of a Riemann sphere and we fix the projection from the double cover to the base by 

: (c + r 1 ) + ^. 0.17) 



4 * i 2 

Conversely, the two pre-images of a point z (jt [a, b] are given by £(z) and £(.z) _1 with 

CW= (^ + ^) a , (3.18) 
o — a 

For all square roots we choose the convention that there is a branch cut from — oo to and that 
VT = 1. Then z ^ [a, 6] implies > 1. In fact, if we write z = x + iy + 2±2 5 the curve 

\C(z)\ = r > 1 is the following ellipse with centre 

x 2 y 2 b—a r 2 + l b—a r 2 — l 

r£ r,i 2 2r 2 2r 

y 

This ellipse can also be described by \z—a\ + \z—b\ = ^^-(r + r _1 ), which shows that a and b are 
the two focal points of the ellipse ((3.19(1 . The significance of this parametrisation is that the contour 
integrals defining 1)1. 1|) will be suppressed by factors r~ 2n when carried out along the ellipse (see 
section l5*T3*|) . 

On the double cover, each J~-field corresponds to the pair of vertex operators V~(()V+((~ 1 ), 
where V±(£) denotes the vertex operator V ±1 ^(() of J 3 -charge ±^g- Thus one expects that 



(n\J (zi)---J (z m )\aa) 



{merer) 



A m 1 

4 ' "' n (^l^-(Ci)^(cr 1 ) • , 



a 



(3.20) 



where Q — C( z i)- The right hand side can be computed in terms of the Coulomb gas expression for 
free boson vertex operators, and one finds 

(nirtzQ-j-^+^Mio) _ / 4 cr 2 " n / C.-C 3 \ 2 ,„ 2n 

(n\a+(b)a-(a)\0) Kb -a) fj & - C 1 ) 2 - ^ ' 

It is straightforward to check that this indeed solves 1(3.1511 : to this end it is useful to rewrite the 
differential operator 1(3.16(1 also in terms of the £- variables, using £ — = rzr-\/ — a)(z — 6) and 
^ = ^^(l — C -2 )^- Of course, 1)3.15(1 only determines 1)3.21(1 up to a constant. This constant can 
be found recursively using 1)3.2(1 . 

(n+l\J-(zi)---J-(z m )\crcr}= [ w 2n+1 (n\J~ (w)J~ (zi) ■ ■ ■ J~ (z m )\aa) pL . (3.22) 

This determines the overall constant for a correlator with m+1 insertions of J~ in terms of a correlator 
with only m insertions. Finally, the expressions with no insertions of J~ has already been given in 
1(3.7(1 . This procedure is the origin of the factor (^) in |3~2Tjt and (|3~2T|) . 

A correlator where some of the J~{zi) insertions have been replaced by J + {zi) insertions can be 
obtained by continuing Z{ through the branch cut [a, b], which amounts to replacing Q — > C" 1 in 

Q3D. 



li 



4 Properties of the operator I\(a, b) 



In the previous section we have collected some information about the structure of the twist field 
correlators. Now we want to study the correlation functions of I\ (a, 6) defined in 

4.1 Monodromy and operator product expansion of sw(2)-currents 

Before computing the monodromy of the sw(2)-currents J , J 3 in the presence of I\(a,b) let us 
consider a slightly more general situation. Define the operator 



B(a,b) = exp^ J F(x)dx 




(4.1) 

where F{x) is a linear combination of holomorphic functions multiplying chiral fields (and we assume 
that the operator product expansion of F(x\) and F(x2) does not have any poles so that no normal 
ordering prescription is necessary). Consider the analytic continuation of a chiral field <j){z) around 
the point b, 

anal, cont. 

a b 



(4.2) 

It is then clear that under the analytic continuation (|4.2|l . the monodromy of 4>{z) is 

<f,( z ) _^ eX p(2« J F(x) ^) 0( z ) . (4.3) 

In the case of I\(a,b) we have F(x) = -^J + (x). Acting on the state |</>) corresponding to the field 
4>(z), the monodromy (|4.3|) then reads \(p} — > exp(iAJ h )|</i). Representing the linear combination 
\<f>) = a\ J + ) + (3\ J 3 ) + 7I J - ) by the vector (a, (3, 7), the monodromy of the sw(2)-currents around 
the endpoint b of I\ (a, b) is given by the matrix 

/l -iX X 2 \ 

M x (b) = Q 1 2iX . (4.4) 

\o 1 ; 

For example, exp(iAJ l ")| J 3 ) = (J 3 j + i\J^ J 3 1 )|0) = |J 3 ) — zA|J + ). Since the currents are single 
valued on C— [a, 6], the monodromy around a is inverse to that around b so that M\(a) = Mx(b)^ 1 = 
M_ A (6). 

The above calculation can be repeated for the stress tensor T(z) and one finds that due to 
[Jq , L m ] = 0, the stress tensor is single valued across [a, b]. Translating this observation back into 
matrix model language (as briefly described in section ll.lf) . gives a way to derive the quadratic loop 
equation of the matrix model from the free boson conformal field theory |2HI 1421 13"3"] . 

To analyse the singularities of the su(2)-currents close to the points a and b it is helpful to 
introduce the following combinations of fields 

J+(z) = J+{z) 

J*{z) = J 3 (z) + Ai nf ^j+( Z ) (45) 

= ^(*)-£k^W-(&k^)V(*). 

Using (|4.4|) one verifies that J ± {z) and J^{z) are single valued in the presence of a single insertion 
of I\(a,b). In particular, they can be expanded in integer modes 

/dz 
(z-p) m r(z)— (4.6) 
. . - 
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Here c S { + , 3, — } and p equals a or 6. We would like to establish the following statement: 

(S) Suppose lim z _>h(2 — b)J~ (z)I\(a, b) is finite inside any correlator, and there exists some N > 
such that lim 2 ^f,(z — b) J + (z)I\{a,b) and lim z _>;,(z — b) N J 3 (z)I\(a,b) are zero. Then 

J^,bh(a, b) = = Jl. b h(a, b) for to >0 , J+^a, 6) = = Z| 6 Z A (a, b) . (4.7) 

Together with the analogous statement for a, this means that the fields at the end points a and b are 
in fact highest weight with respect to the J c action. 

The vanishing of the J + -zero mode implies in particular that J + (z)I\(a,b) is regular for z — ► b. 
Of course, this also follows immediately when writing out the integrals in (|1 . 1|) and using that the 
operator product expansion of J + with itself is regular. However, in section IK. 51 we will need to apply 
statement (S) to S\(a, b) instead of I\ (a, 6), so we will present the argument in a form which will be 
valid also then. 

To establish (S), start by expressing J~ in terms of the single valued combinations l|4.5[) . 

J-(z) = J~(z) + |A In |=| J3(z) - (A i n ^fj+{z) . (4.8) 

Denote by C(<j)) a correlator involving <f> and I\ (a, 6), as well as any number of other fields. Expanding 
the right hand side of (|4.8|l in terms of modes around b then gives 

N 2 

C(J-{z))= ^-Q^faj-J + ghL&C^) - (&\n&) <?(J+ fe )). (4.9) 

m— — oo 

The summation is truncated by the assumption on the limit z — > b of J + and J 3 . Evaluating the 
conditions lim z _^£,(z — b) m+1 J~ (z)I\(a, b) = for m = N, N—l, . . . , 1 gives C(J _;?) = for that range 

of to. Finally, for to = we get C( J^ 6 ) = = C(J$. b ). Since C(. . . ) was an arbitrary correlator, this 
implies (|4.7j) . 

In order to establish the highest weight relations 14.7fl for I\(a, b) we still have to verify that the 
conditions of the statement (S) are met. For J + this is obvious, but for J 3 and J~ this requires a 
short calculation which is given in appendix lA.il 

Finally, let us show that the endpoints of I\(a, b) obey the Virasoro highest weight condition for 
weight zero. To this end, instead of expressing the stress tensor T(z) as in Ij2.3|l we use the single 
valued fields (|4.5I) . Computing the operator product expansion of J 3 with itself to order 0(z— w) one 
finds 

T(w) = lim (J^JV) - 7^^) + & 7" ,J + (w) , (4.10) 

z^w V (z— w) z / z7r [w — ajyw — b) 

which in terms of modes around b reads 

A °° 

L m , h = ■ Jl- b Jl-k-b ■ +^ 5](a-&r fe J+ +k ; b ■ (4.11) 
fcez k=o 

Together with (j4.7|l this immediately implies that L m -i,I\(a, b) = for to > 0. For a instead of b one 
finds the same result. 

4.2 Zero-point function 

Up to now we have kept the parameter A arbitrary. But just as was the case for twist fields, we can 
restrict our attention to the case A = 1. The results for general values of A are then obtained via the 
identity 

h (a, 6)|0) = e ln W J o Zito, 6)|0). (4.12) 

To see this write I\(a,b) and Ii(a,b) as a sum over J + integrations and use [Jg, J + (z)] = J + (z) 
which leads to exp(ln(A) Jq ) J + {z) = XJ + (z) exp(ln(A) Jq ). Moving the exponential in (|4.12|l past the 
J + insertions in the expansion of I\ (a, b) gives a factor of A for each insertion. 
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The correlator we would like to compute is (n\Ii(a, b)\0). From the previous section we know that 
the endpoints of I\(a,b) behave as Virasoro primary fields of weight zero. Applying a rescaling and 
a translation yields 

(n\h (a, b)\0) = {b-a) n2 2- n2 (n\h (-1, 1)|0> . (4.13) 
The correlator on the right hand side can be computed by explicit integration using orthogonal 
polynomials on the interval [—1, 1], i.e. the Legendre polynomials (see e.g. ^J|22] for an introduction 
to the method of orthogonal polynomials), or by directly using Selberg's integral (see e.g. |371 chapter 
17]). The result is 



1 

(2tt)" n\ 
1 



<n|Ji(-l,l)|0) = j dx n (n\J+(xi) ■ ■ ■ J+(x»)|0) 



(4.14) 

J dxx--- dx n A(x) 2 = (2tt)"" 2™ 2 det(H n ) . 



(27T)" 

In the first step, the definition l|l.l|l has been substituted, in the second step the Coulomb gas 
expression for the integrand has been written in terms of the Vandermonde determinant A(x) = 
Y[i>j( x i~ x j)- I n t ne result, H n is the nxn-Hilbert matrix (H n )ij = (i + j — 1) . Its determinant 
is given by 

, , v 2 '-r-j 2 (2 k k\ 2 \ 2 2"- 2 « 2 7 r«+ 1 /2 / G ( n+ i) G (i)\ 2 

Here G(z) is the Barnes function, which is defined by G(z+1) = Y{z)G{z) and G(l) = 1 together 
with a convexity condition. The behaviour of det(H n ) for large n can now be obtained from the large 
z expansion of G(z). The latter can be found in ^ eqn. (28)] (in the preprint (vl), the minus in front 
of the sum should be a plus) or in [251 ecm - (2.38)]. In this way, we finally obtain 

{n\Ix(a, b)\0) = (6-a) n V n V4 2 A e 3C'(-i) exp ^ _ i_ n -2 + _i_ n -4 + ( n -6^ _ ( 416 ) 

Comparing to the correlator of two twist fields l|3.7|l we observe that a + (b)<j-(a) does correctly 
reproduce the leading term in l|4.16(l in the sense that 

lim ^ 2 ln(n|cr + (&)cr_(a)|0) = lim vT 2 ln{n\h(a, 6)|0) . (4.17) 



4.3 One-point function 

After computing the zero-point function we will turn to the correlators (n\ J c (z)Ii(a, b)\0) for c € 
{+, 3, — }. The result can again be obtained by matrix model techniques using orthogonal polynomials 
(see e.g. chapter 22] or |221 section 4]), which amount to explicitly computing the relevant 
integrals. Another method more in the spirit of conformal field theory is to use that the h = 1, c = 1 
Virasoro highest weight representation has a null vector at level three. Denoting the highest weight 
vector in this representation by | J) the null vector \rj) is 

\rf) = (3L_ 3 - 2L_ 1 L_ 2 + J) - . (4.18) 

Since the three su(2) currents J c (z) are Virasoro primaries of weight one, the three correlators 
(n| J c (z)/i(a, b)\0) will all solve the same third order differential equation in z, obtained from the 
null vector \rj). To compute this differential equation, recall from section [3. II that from the point 
of view of the Virasoro algebra there is no difference between an insertion of I\ (a, b) and a product 
<p(b)<t>(a) of Virasoro primary fields <fi with conformal weight zero. Specialising to a = —1 and 6=1, 
the differential equation is then found to be (see |121 section 8.3] for more details on null-vector 
computations) 

= <n|»j(2r) #1)^-1)10) 

f 1 4z 5z 2 -l 2(n 2 -l) / z \ 1 ( 4 - 19 ) 
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The space of solutions to this equation is three-dimensional. The elements in this space describing 
the three functions (n| J c (z)Ii(— 1, 1)|0) have to be identified from their behaviour at the singular 
points —1, 1, oo. 

In any case, using either orthogonal polynomials or Virasoro null vectors, the final result for the 
one-point functions is 

(n\J+(z)h(-l, 1)|0) = 7rn[Pn-i(z)P^z) - P n {z)P' n _ l {z)) (n\h(-l, 1)|0) 

(n\J s (z)h (-1,1)|0) = n(p n _ 1 (z)Q;( 2f )-P n (z)Q;_ 1 («))(n|J 1 (-l,l)|0) (4.20) 

(n\J-(z)h(-l, 1)|0) = - -(Qn-xWQ'^z) - Q n (z)Q' n _ 1 (z)) (n\h(-l, 1)]0> . 

In these equations, P n (z) is the n'th Legendre polynomial and Q n (z) the n'th Legendre function of 
the second kind. The first few are (chosen such that they are real for z £ M. — [— 1, 1]) 

1 / z+1 \ , . z , / z+l \ , , 3z 2 — 1 / z+1 \ 3z 

One can verify that the functions (|4.2U|) solve l|4.19|l and their monodromy around the point 1 is given 
by (|4~4l . 

To understand the large- n behaviour of the one-point functions, it is convenient to write P n (z) 
and Q n {z) for z £ [— 1, 1] in terms of hypergeometric functions as [23 eqn. (8.723)] 



I>+|) C" +1 

?r(n+i) Tc^T 31 ' u -'^ : '~ 



P„(z) = c ^^ -Wt^' 1,1 



■2> ]\ _ 

(4.22) 



0Fr(n+i) c- wi 1>3 . _x \ 
-* ~ r ( n+ 3) y^fzf 2 l( ^2' 2' 2+ n > c^tJ ■ 

Here, C = C(^) is given by l(3.18|l with a = — 1, b = 1. The large- n expansion of the correlators 1(4. 20|) 
is then found by writing out the definition of the hypergeometric functions ((4.22(1 as a power series 
and expanding the Gamma- functions, 

(n\J+{z)h (-1,1)10) _ 2C 2 " ( l C + r 1 -x . 1 (1+5C 2 )(1+5C- 2 ) 2 , 3 



(n|J x (-1,1)|0) (C-C- 1 ) 2 V 4 C-C- 1 32 (C - C" 1 ) 2 

(nlJ 3 ^)/^-!,!)^) _ 2n / 1 n 

H/r (-1,1)|0) " C-C- 1 ^ 2CC-C- 1 ) 2 1 ^ 

(n\J~ (z)h{—l, 1) [0) 2r 2 " / IC + r 1 _! , 1 (l+5C 2 )(l+5r 2 ) _ 5 



HJx(-l,l)|0) (C-C- 1 ) 2 ^ 1 4C-C" 1 " 1+ 32 ' "(C-C"' 1 ) 2 ' ' n ~' + ("~ d ), 

(4.23) 

Comparing to l|3.9|l and (|3.21|) (specialised to a = —1 and b = 1) we see again that the leading 
behaviour in 1/n is the same for Ii(a,b) and the twist fields er + (fr)cr_ (a). 

In fact not only the leading term in the expansions (|4.23|l has the Z2-symmetry J ± — > J T and 
J 3 — > — J 3 upon analytically continuing z through the branch cut [—1,1], but this monodromy is re- 
tained at any finite order in the expansions (|4.23|l . It is only after summing all terms that the 'correct' 
monodromy 1(4.4(1 is recovered. For example, at finite values of n, the correlator (n| J + (z)I\{— 1, 1)|0) 
is a single valued function of z (in fact, a polynomial), while in the 1/n-expansion it has a branch cut 
on [—1,1]. Recall from 1)1.7(1 that this conformal field theory correlator is related to the correlator 
Z% e U\det(z - Af) 2 ](") in the hermitian one-matrix model, and hence the latter shares with the corre- 
lator of J (z) the property that the monodromy of the individual terms in the 1/n-expansion differs 
from the monodromy of the complete expression. This is a manifestation of Stokes' phenomenon as 
mentioned at the end of section fTTTI 



5 Writing I\(a,b) in terms of twist fields 

In sections[3]and0]we have collected some properties of sw(2)-twist fields and of the operator I\(a, b). 
We have seen that in sectors of large J 3 -charge, I\{a, b) behaves very similar to a product <7 + (5)cr_ (a) 
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of twist fields. Motivated by these observations one can now try to find an alternative expression for 
the operator I\(a, b) in terms of appropriately dressed twist fields. 

To achieve this, we make an ansatz S\(a, b) involving twist fields and ./"-integrals in section 15. II 
and show that it can reproduce the correct monodromy. That the su(2)-currents have the same 
monodromy for I\ (a, 6) and S\(a. b) is the first piece of evidence for the proposed identity I\ (a, b) = 
S\(a, b). The J~ integrals in S\(a, b) are divergent and need to be regulated; this is done in section l5~21 
Next, in section T5. 31 we investigate the large-n behaviour of S\{a, b) and find that it is given in terms 
of a product of twist fields <7 + (5)cr_ (a). This is the second piece of evidence for I\(a,b) = S\(a,b), 
since we saw, in sections 14.21 and 14.31 that the large-n limit of the zero- and one-point function 
in the presence of I\(a, b) agrees with the corresponding twist field correlators. The requirement 
that S\(a,b) should have the correct monodromy still left an (a, 6)-dependent normalisation factor 
undetermined, which is fixed in section f5. 41 by matching it against the leading large-n behaviour of 
(n\I\(a, b)\0). Finally, in section f5.5l we present the third supporting evidence by checking, to the 
extend that we were able to calculate it, that the su(2)-currents have the same singularities in the 
presence of S\(a, b) as were seen for I\ (a, b) in section l4~D 



5.1 Reproducing the monodromy of I\(a,b) 

Our first task will be to modify the product a + (b)(j_(a) in such a way that instead of the Z 2 - 
monodromy, we find the monodromy 14. 4|) for the su(2)-currents. As a motivation for the ansatz 
below, compare the definition (|l.l|l of I\ (a. b) to the expression (|2.20() for the product a + (b)<j-(a). It 
seems one can go from the latter to the former by 'subtracting' the J~ -contribution to the integrals. 
In this spirit, define an operator 

S(a,b) — C(a,b) a+ 7 (b)exp( — f J^(x)dx S ) exp (— [ J~(x)dx] cr_ 7 (a) . (5-1) 

L V 27T J Cl / V 27T J C2 I J rcg 

The contours Cyp are as shown in 1(1.12(1 . C(a,b) is a C-valued function and a, f3, 7 are constants. 
Let us postpone the discussion of the regulator to section 15.21 in any case it will be defined in such a 
way that it does not affect the monodromy. 

The monodromy of the su(2)-currents around the point b of S(a,b) is a product of three terms, 
two from the exponentiated ./"-integrals and one from the Z 2 -monodromy of the twist fields. In 
terms of the basis used in (|4.4|l the monodromy is 

/l 0\/0 7 2 \/ 1 0\ / aV iaj 2 7 2 

M{b)= -2i/3 1 -10 -2ia 1 = 2ia(l-af3j 2 ) 2a/3 7 2 - 1 -2i(3j 2 
\[3 2 i/3 l/\7" 2 /Va 2 ia l) \ 7 - 2 (l-a/?7 2 ) 2 i/3(a/3 7 2 -l) /3 2 7 2 

(5.2) 

We see that this matches with the monodromy M\(b) of I\ (a, b) in l|4.4|l for precisely two choices of 
the three parameters a, f3, 7, namely a = (3 = —A" 1 and 7 = ±A. We will choose 7 = A. 



5.2 Regulating the J -integrals 

To regulate the J" integrals in 1)5.1(1 we will first introduce a cutoff e and then present a subtraction 
scheme which we conjecture to give a finite e — > limit. In fact, we will regulate ((5.1(1 for the choice 
7 = 1, which will give the operator S\(a,b) in 1(1.11(1 for A = 1. General values of A will then be 
obtained as in 1(4.12(1 . 

The e-cutoff is imposed simply by changing the integration contours to approach the points 
a, b only up to a distance e. More precisely, fix a small positive constant A as well as a value e A. 
Consider the integration contours C\ t i and C e a h defined as follows, 




1G 



The dashed circles around the points a and b have radius e so that C E a — [a— A, a— e] and C[ = 
[b+e, fr+A], with orientations as indicated. Indeed, for e = this is just a deformation of the contours 
Ci,2 defined in l|1.12(l . Instead of integrating J~(x) define a field pt{x) as 

Pt (x) = tJ-( X ) - m{ J-(x))l , where . ^'^V^t^g 1 ^ ■ (5-4) 

(0|<t + (6)ct_(o)|0) 

Here t is a complex parameter and /(i) = t + f2t 2 + f^fi + . . . is a power series in t. On the segments 
and C1.2 define the following integrated fields, 

U+j{b) = exp(-^- J^p t (x)dxja + (b) 

U^j(a) =exp^-iy p t {x)dx^a^{a) (5.5) 
v t A (. a > b ) = C( a > b) exp ( - ^ pt(x)dxj exp ( - ^- ^ />t(x)da;) . 

In there is a factor of i instead of because in deforming the contour (|1.12() to (|5.3|) . the 

segments b are traversed twice. 

The subtraction scheme now consists of expanding the operators U± t as a formal power series in 
t and demanding that each term in the expansion has a finite limit as e — > 0. This procedure will 
result in conditions determining the constants /a, /3, .... It is not at all obvious that such a solution 
exists, and we have no proof that it can be done to all orders in t. In appcndix lA.2l we verify that the 
subtraction scheme ()5.5JI works at least to order < 3 , with f(t) = t — ^t 2 + i< 3 + 0(t A ). To proceed 
we will assume: 

(Al) There exists a function f(t) such that each order in the expansion of the operators U+^(b) and 
Utl\{a) in powers of t has a finite limit as e — * 0. 

Using (Al), we can define the operator S%(a, b) in terms of U±\ as 

Sx(a, b) = Vf{a, b) S p + e £(b) S p _ e £(a) , (5.6) 

where 

Sl e X(b) = (limt^>)) t=1 , STfca) = (limC/^(a)) t=i . (5.7) 
Finally we then obtain 

S x (a 1 b)\0)=e h ^S 1 (aM0)- (5.8) 

A number of comments are in order. First, the abbreviation 'pert' in l|5.7|l stands for 'perturbative', 
a qualifier that will be justified in the next section. Second, while the decomposition l|5.6|l of Si(a, b) 
will be useful in the following, we can equivalently write it in a form that more closely resembles 

Si(a.b) = C(a, b) ^ lim a + {b) exp ( - ^- J ^ p t (x)dxj exp ( - ^- J p t (x)dx^a^{a) S j . (5.9) 

Here, Cf is the contour obtained by joining C^, C\ and Cf , and Cf is obtained by joining C^, C2 and 
CJ. In the form l|5.9|l it is also apparent that S%(a,b) has the monodromy l|5.2[l . since pi(x) differs 
from J~ (x) only by a central term. 

Third, it is clear from the definition that S^^^b) and 5^ e ^(a) are coherent states in the Z2 twisted 
representations generated by a + and er_ , respectively. Coherent states in twisted representations also 
appeared in relation to matrix models in . 
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5.3 Behaviour of S\(a,b) at large n 

In section f4.3l we have seen that the monodromy of I\(a,b) around a or & is entirely due to 'non- 
perturbative' effects, i.e. that to any finite order in l/n, the monodromy is just given by the Z2 twist 
H1.9(l . We now want to show that the behaviour of S\(a, b) is the same; more explicitly, we shall show 
that (see QTTty ) 

(n|(ficlds)5 A (a,6)|0) = (n|(fields)Sf unc (a, &)|0) (l + 0(r~ 2n )) , (5.10) 

for some r > 0, where S\f unc (a, b) has the same monodromy as the two Z2-twist fields. We will again 
only treat the case A = 1; for the general case the reasoning is analogous due to (|5.8|l . 

To define S\ rnnc (a,b) we choose the contours C\ : 2 in i|5.3|) to lie on the ellipse (|3.19() of constant 
|£| passing through a — A and b + A. This amounts to choosing r = |£| to be 

r = l + a?(l+ V 1 + ^)- ( 5A1 ) 

Then define the function 

D A (a,b) =C(a,b)exp^f(t) J_ (J~ (x)) dxj , (5.12) 

and set 

5i runc (a, b) = D A (a, b) S*$(b) S^(a) . (5.13) 

The claim H5.1()|l is then implied by the following statement: let (F) = Ill=i J Ci ( w i) abbreviate a 
product of sw(2)-currents. Then we have 

(n\{F)(V A {a,b)~ Dt(a,b)l)r)+(b)r)-(a)\0) = (n\(F)r] + (b)r)-(a)\0) ■ 0(r~ 2 ™) , (5.14) 

where r]±(z) are fields corresponding to states in the Z2-twisted representations. This is sufficient to 
establish l|5.10|l since the fields S^ Ji (z) correspond to (coherent) states in the Z2-twisted representa- 
tions (and thus can play the roles of r)±). 

To prove l|5.14|) . first note that it is enough to consider products (F) = Y[[=\ J Ui (wi) with i/j = ±, 
because J 3 insertions can be obtained in the operator product expansion of J + and J~ . To express 
twist-field descendents, define the modes M rjS as 

M rs = f (z-a) r (z-b) s J a (z)—, r,seZ+± (5.15) 
Jc ab 2tti 

where C a b is a contour encircling a and b. Since we are at level k = 1, the entire Z2-twisted repre- 
sentation is generated by acting with modes J r 3 , r £ Z<o+5 on cr±. Correspondingly, the product 
?7+(6)?7_(a) can be written as a linear combination of appropriate products fL M riSi a + (a)a-(b), with 
fi+Si < 0. Next, note that by definition, 

00 1 t m f 

V A {a,b)- D A (a,b)l = D A (a,b) V — (^-) / J~(xi) ■ ■ ■ J~{x m ) dxi ■ ■ ■ dx m . (5.16) 

Upon inserting (15 . 1 ti|) into 1)5. 14[) one obtains integrands of the form 

M 

I m = (n\(F)J-(xi)---J-(x m )Y[M ri , Si \aa), \aa) = o+{b)o-{a)\0) . (5.17) 

i=i 

Using the commutator [M rjS , J ± (a;)] = ±(x— a) r (x~ b) s J* (x), as well as (n\M r ^ s = S r+S fin(n\ (expand 
M r>s in integer modes of J 3 and use the condition r+s < 0) we can write 

I m =A(n\(F)J-(x 1 ) ■ ■ ■ J-(x m )\aa) 

M \F\ m (5.18) 

A = Y[ ( - v 3 {w 3 -a) r ^{w 3 -b)^ + J2( x j- a ) n ( x j- b ) Si + nS ri+Si fi) ■ 
*=i 3=1 3=1 
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To proceed we need the following two estimates on the large n behaviour of correlators, 

= (const) + 0(n~i) , = . 

(5.19) 

For the first estimate note that the correlator in the denominator produces a factor given by A as 
in H5.18fl but with m = 0. To see the second estimate, insert the explicit solution <|3.21[1 for the 
correlators. For each J~(xk) insertion there is a factor of ((xk)~ 2n in the numerator which, since the 
Xk lie on the contour Ck, has = r. Then 

T tm\(w\ TT M I X A(n\(F)\aa) (n\(F )J~( Xl ) ■ ■ ■ J-(x m )\<ra) 

J m = (n \(F) I I M ruSi \aa) rrE^TTTl i — 7 7 \ 

11 (n\(F)U M r„sMa) (n\{F)\cra) ( 5 20 ) 

= (n\(F)l[M ri , Si \aa) ■ 0(r- 2nm ) . 

When inserting (|5.16l) into l|5.14[l . the most relevant contribution therefore comes from m = 1 and 
thus is of order 0{r~ 2n ). 

5.3.1 Approximating S\(a,b) in terms of twist fields 

Next we want to show that to leading order in 1/n, S\(a,b) can be replaced by a product of twist 
fields (see l|1.17|l ). To obtain this relation we start by defining a field pt, n {x) in the same way as 1)5.4(1 . 
but where we subtract the one-point function with respect to (n\ instead of (0|, 

Pt , n (x)=tJ-(x)-f(t)(J-(x)) n l, where (J-(x)) n ^ i^^^figM. , (5 .21) 

Let C be any contour from a to b. The difference pt{x) — pt. n {x) has a well-defined integral along C 
since the singular contribution from the poles at a and b cancel. Using 1(3.21)1 one finds explicitly 

r r / — 2n i ri ~ 1 o 

/ (p t (x)-p t , n (x))dx = f(t) ^ d( = f{t) y—7 

= /(t)(V(n+|)+7+21n(2)) ( 5 - 22 ) 

- /(t)(ln(n) + 21n(2) + 7 + ^ + 0(n^ 4 )) , 

where C is any contour from —1 to 1 not passing through £ = 0, 7 is Euler's constant and ip(z) is 
the digamma function, ip{z) — T'(z)/T(z). We can then consider the product 

O n (a,b) = exp(^- [ (ptM-pi^fidxjSifab), (5.23) 
KZlT Jd+c 2 ' 

which amounts to replacing p\(x) in the definition (|5.9|l of Si(a, b) by pi iTl (x). When using p f rl instead 
of pt, both the ./"-integrals and the integrals of the one-point functions (J~(x)) n are suppressed away 
from a and b for large n. It is then plausible that in the expansion of the exponential (|5.9|l . all terms 
involving /9t,„-integrals are of order 0(n _1 ). This can be checked explicitly for the first few terms in 
the expansion, but we have no general proof. Let us hence assume 

(A2) The operator O n (a, b) in l|5.23|l has the large n behaviour 

(n\ (fields)C„ (a, b) |0) = C(a, b) (n\ (fields)<7 + (b)a- (a) |0) (l + 0^ )) . (5.24) 

In order to obtain a 1/n-expansion, one should thus not consider the correlator (n\ (fields) S\{a, b) |0) 
directly, but rather the normalised version (n|(ficlds)0„(a, 6) |0)/(n| (fields) |ctit). However, it is not 
true that the expansion of the exponential l|5.9|l with p t ^ n instead of p t produces the l/n-expansion 
term by term. Instead even a term with m J~ -integrations gives a contribution of order n . This is 
not so surprising when one considers more carefully the regulated expression for the m'th term in the 
expansion. In fact, the coefficient of t m will also contain a term of the form f m J (J~ {x)) n dx from 
the ^-expansion of J p t , n (x)dx. In this sense the subtraction scheme mixes all orders and it is only 
easy to extract the large-n limit, but not the subleading terms in the l/n-expansion. 
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5.4 Fixing the normalisation of S\(a,b) 

In section lh. II we have partly fixed the operator Si(a, b) by requiring it to have the same monodromy 
as 7i(a, b). Using assumption (A2) we can further fix the normalisation C(a,b) by demanding 
(n\Si(a,b)\0) / (n\li(a,b)\0) = 1 in the largc-n limit. To this end, combine ((5. 22(1 and (|5.23|) to 
obtain 

O n (a, b) = n-e^WMS^a, b)(l + Ofa" 1 )) . (5.25) 
Next, using the above approximation together with 15.24f) and ((3.7(1 leads to 

{n\Si(a, b)\Q) = C(a, 6)n-^ i e -^ i(21n(2)+7) 4-" 2 (6 - a)" 2 -* (l + O^ 1 )) . (5.26) 
Comparing to the corresponding expression 1(4.16(1 for l\ (a, b) then leads to the requirements 

/(l) = f , C(a, b) = 2T2 e 3C'(-i)+7/4 (6 _ a) | _ (5 27) 

For course, f(l) is in principle determined by demanding the e-limit in (|5.9(l to exist, but one would 
need the expansion to all orders to check whether its values is indeed f. It may be taken as an 
encouraging sign that from f(t) = t — ^t 2 + ^t 3 + 0(t 4 ), the first three approximations to f(l) — j 
are 0.21, -0.10, and 0.06. 

After fixing C(a, 6), the operator Si (a, b) is completely determined. The definition of the regulated 
expression (f 1 . 1 1|1 is that Si(a,b) is given by 1)5.9(1 with C(a,b) as in ((5.27(1 . In 1(1.11(1 . the constant 
multiplying (6— a) s in C(a, b) has been absorbed into the definition of [. . . ] rcg . 



5.5 Singularity structure of the s-u(2)-currents in the presence of S\(a,b) 

Finally we want to argue that the leading singularities as the currents approach the endpoints of 
Si(a,b) are the same as for the case of I\{a 1 b). For 7i(a, b) these singularities were analysed in 
section l4~Tl and the relevant properties are summarised in the statement (S). Here we want to present 
two pieces of evidence that the conditions of the statement (S) are also met for Si (a, 6) in place of 
Ji(a, b). In particular we want to argue that the leading singularity as J~ (z) approaches the endpoints 
of 6*1(0, b) is a simple pole. We start by investigating the behaviour of Si (a, b) under global conformal 
transformations . 

Let (p(z) = ^fqi^ be a Mobius transformation and U v be the operator implementing that trans- 
formation on the space of states (an explicit expression in terms of Virasoro generators can for 
example be found in (241 section 3.1]). Denote the field pt(z) introduced in ((5.4(1 by p t (z;a,b) to 
keep track of the values for a and b entering its definition. Using U V J~ (z)U v -i — (p'(z) J~ (ip{z)) and 
U ip a±(z)U v -i — ip' {z)ts a±(ip(z)) it is easy to check that 

U v p t (z; a, b)U v -i = ip'(z) p t (ip(z); ip(a), ip(b)) . (5.28) 

Next we observe that for C(a,b) defined in 1(5.27(1 . we have C(a,b)(ip' (a)<p' (b))™ = C (ip(a) , tp(b)) . 
Thus the total effect of the Mobius transformation on Si (a, 6) is 

U v Si{a,b)U v -i = C(a, b) I lim U v exp ( - — f p t (x; a, b)dx)a + (b)a-(a)U ip -i ] 

= Si(^a),^(6)). (5.29) 

We see that Si (a, b) transforms as a product 4>{a)4>(b) of two Virasoro quasi-primary fields of conformal 
weight zero. This observation allows us to compute the correlator (0| J c (z)Si(a, b)\0), for c e {+, 3, — }, 
since conformal invariance fixes a three-point function up to a constant, 

(0| J c (z)S 1 (a, b)\0) = (const) (z-a)- 1 {z-by 1 {a~b) . (5.30) 

In particular this shows that for the out-state (0|, all su(2)-currents have at most a first order pole 
at a and b. 
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As a second piece of evidence, we compute the correlator (n\ J v (z)S%(a, b)\0) for v = ± to first 
order in the t-expansion (that is, in the definition l|5.9[l of Si (a, b) we expand the exponential to first 
order in t before setting t = 1). Instead of using Si (a, b) it is convenient to use O n {a, b) as defined in 
H5.23|l which differs from Si(a,b) by a constant. One finds, with £ = C(z) and £ = ((x), 

(n\r(z)O n (a,b)\Q) 

= C(a,b)({n\r(z)(l-±. [ {j-(x)-(J-(x)) n )dx + 0(t 2 ))\aa)) (531) 
\ JCi+c 2 j t=1 

= C(a, b){n\r{z)\aa) (l + 1) jf ^ ( ^f!" 1)2 ^ + °^)) ■ 

Here C[ is a contour from —1 to 1 passing along the upper half of the unit circle, while C' 2 passes 
along the lower half of the unit circle. The integral in (|5.31|) is given by 



■ dx = 



c> {xy — l) 2 L 2n + 1 (xy — 1)' 



2 F 1 (l,2;2n + 2;(l-xy)- i ) , (5.32) 



x=l 



x=-l 



where C is a contour from —1 to 1, which also determines the relevant branch of the hypergeometric 
function. We are interested in the asymptotics of l|5.31(l for ( — > ±1. This amounts to taking the 
argument of the hypergeometric function in H5.32JI to infinity. The hypergeometric function has the 
asymptotics, for u — > 00, 

2^(1, 2; 2n+2; u) = -(2n+l) u" 1 + 0(ir 2 ln(u)) . (5.33) 

Altogether we find that the leading singularities of the integral (|5.32l) are first order poles at y — ±1, 
which in l|5.31|l get cancelled by the prefactor C^ v — 1. Hence to first order in t, the leading singularity 
of (|5.31(l for z — > a, b is that of (n\ J u (z)\aa), which clearly satisfies the conditions of (S). 

After presenting these two calculations regarding the poles of J c {z)S\{a, b) for z — > a,b, we will 
assume that in general 

(A3) for x — a, fe, inside any correlator, lim z ^ x (z — x)J~ (z)S\(a, b) is finite, and there exists some 
N > such that lim z _^ I (z - x) N J + ' 3 (z)Si(a, b) is zero. 

It then follows that the fields at the endpoints of Si (a, b), just as for ii(a, b), are highest weight with 
respect to the single valued combinations l|4.5|l . and that in particular they are Virasoro primary of 
weight zero, i.e. L m . a S\(a, b) = = L m -b S\(a, b) for m > 0. Together with the fact that I\ (a, b) and 
S\(a,b) have the same monodromy properties, this is very good evidence for the equality of I\(a,b) 
and S\(a, b). 



6 Outlook 

In this paper we have proposed and provided evidence for the operator identity I\{a,b) — S\(a,b). 
Here I\(a,b) has a simple formulation as an exponentiated integral of J + -currents and is directly 
related to eigenvalue integrals in matrix models. The expression for S\(a,b) is more complicated, 
involving twist fields and a regulator. However, properties of the large-n limit of correlators are easily 
understood in terms of S\(a, b) while they are harder to see when using I\ (a, b). 

There are several directions in which one can attempt to generalise the analysis of this paper. 

(i) From the conformal field theory point of view, the most obvious question is whether there is a 
generalisation to su(2) at level k > 1. In this case the operator I\(a, b) is still given by Ijl.lfl . However, 
there are now fc+1 Z2-twisted highest weight representations, and there is no longer a unique (up to 
scalar multiples) out-state (n\ which is highest weight with respect to and has J 3 -charge n, but a 
finite-dimensional subspace. Both effects make it more difficult to identify the analogue of S\(a,b). 
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Independently, it is also of interest to see if correlators of I\(a, b) in su(2)k do have a matrix model 
interpretation. 

(ii) Recently a conformal field theory approach to the non-linear a model with the complex torus 
(C*) d as target space has been investigated in It uses a non-unitary CFT which contains a su(2) 
algebra at level k — 0. Similar to (|f .ljl it has an integrated exponential operator, which does not 
require normal ordering and generates a logarithmic branch cut. The above model is related to the A- 
model by a deformation and to the B-model by an additional T-duality. It would be very interesting 
to understand whether this CFT approach to the non-linear a-model is linked to a matrix model 
description. On non-compact Calabi-Yau 3-folds associated to ADE singularities [13] the topological 
B-model, a subsector of the non-linear cr-model, has already been related to ADE-quiver matrix 
models 11411151 17). see also (iv) below. 

(iii) From the matrix model point of view, one should use more general potentials than the infinite well 
potential that we considered here. This amounts to taking the out-state (n\e~ H instead of just (n|. 
Also, it would be good to remove the effect of the hard edges, possibly by considering S\(a — e n , b+e n ) 
to shift the endpoints by an n-dependent amount away from the location of a cut. 

(iv) There is a whole class of multi-matrix models, called ADE-quiver matrix models, which can 
be rewritten in terms of free bosons, or, more precisely, in terms of an ADE-WZW model at level 
one [311 1291 1421 1321 1151 [7J. For these one would define a number of operators P x (a,b), indexed 
by simple roots on, i = 1, . . . ,r. Since each of the corresponding raising operators E ai lies in an 
sw(2)-subalgebra, one would expect that the analysis of this paper can be repeated without much 
modification. The comparison of perturbative versus non-perturbative moduli spaces of FZZT-brancs 
(as mentioned in the introduction) has also been carried out for (j>, l)-minimal string theory using a 
two-matrix model |28| (however, not an ADE-quiver model) and it would be interesting to compare 
results. 

(v) Insertions of several operators /^(ai, £>i), I\ 2 (a,2, 62), • ■ ■ correspond to several cuts in the complex 
plane. On the matrix model side one obtains in this way not a partition function with fixed filling 
fractions, but rather a generating function in the parameters for the various filling fractions. It 
would be interesting to see if the methods of this paper can be extended to be a useful tool in 
investigating the 1/n expansion of such multi-cut solutions. 

We hope to address some of these points in the future. 
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A Appendix 

A.l Calculation of J 3 and J pole with I\(a,b) 

In this appendix we show that the J 3 - and J~ -conditions of statement (S) in section f4.1l are met. 
For the case of J~ we need to prove that 

lim (z - 6) ( TT J ai (wj) J-(z)I x (a,b)) = finite, (A.l) 

z^b * XX 

3 
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where the Wj ^ a, 6 are pairwise disjoint. Then 
(l[j a >(v,j)J-{z)I\(a,b)) = ±- l J^dz x ■■■ J^dz n (Y[j^(w J )J-(z)J + (z 1 )---J+(z n )) 

_ (A.2) 

for some suitable n. The amplitude in the integrand can be calculated recursively, using the singular 
part of the operator product expansion (see for example Starting with the J + fields, it is then 

obvious that the terms that could be singular in the limit z — ► b arise in two ways: either we have 
the double pole 

J + ( Zi )J-(z)~ 1 (A.3) 
\z% z J 

that gives rise, after integration of Zi, to a simple pole in (z — b); this contribution is therefore finite 
in the limit of IjA.ljl . The second contribution conies from the simple pole 

J+( Zi )J-(z)~^M. (A.4) 

Zi~ z 

The other J + (zj) fields can then either contract with the J aj (wj) fields, or we can get a further 
contraction of J + (zj) with J 3 {z), which then leads to 

J+iz^J+iz^J-iz) ~ - (A.5) 
(z z%)\z Zj) 

In either case, it is straightforward to determine the integrals, and we obtain either a log(z — b) or 
a log 2 (z — b) term. In the limit of i|A.l|) these contributions therefore vanish. 

The analysis for the case of J 3 is essentially identical; in this case, only the second type of terms 
appear, and we find with the same arguments as above that 

lim - b) ( TT 3 a i (wj) J 3 {z)I x (a, b))=0. (A.6) 

3 

A.2 Existence of the regulator to order t 3 

In this section we compute the first few orders in t of the function f(t) = fit + f^t 2 + ht 3 + 0(t ) 
which enters the definition l|5.6|> of Si(a, b) via l|5.4|l . To this end we make use of the decomposition 
H5.6|) of Si (a, 6). The operators £±°a are well defined if U± t (as given in 15.5(1 ) has a finite e — > 
limit. This requirement fixes the constants /i, /2, /3 uniquely. Here we will only treat 5?^ (6) 
explicitly. The calculation for S^ e A (a) is analogous and leads to the same answer. 

In section f5.5l it was shown that under global conformal transformations, S\(a,b) behaves as a 
product <j)(a)4>{b) of Virasoro-primary fields of weight zero. We will use this freedom to assume that 
6 = and a = — oo. The question whether S^^(b) is well-defined now amounts to checking that 

C/?£(0)|0) = cxp (i J p t {x) dx) \a + ) (A.7) 

has a finite e — > limit, order by order in t (the sign difference with respect to (|5.5H is due to the 
change of direction in the integral). Here, pt{x) takes the form 

_ » , . (0\J-(x)a + (0)a-(a)) , x fit) lk . 

Pt (x)=tJ-(x)-f(t) lim V ' , V ' \ " = tJ-{x)-^. A.8 

a^-oa (Q\<j + (V)<j-{a)) 4x 

Set further R(x) = J~(x) — fi/(Ax). Then the first few orders in the i-expansion of 1A.7|I read 
U E + ^(0)\0)= \<r+)+t- I R{x)dx\a + ) +t 2 (\ f R(x)R(y)dxdy-- / "A dx)\a+) 

TT J VtH J x>y TT J AX J 



+ t 3 (^f R(x)R(y)R(z) dx dy dz - \ f R(x)j-dxdy (A.!)) 

n Jx>y>z n Jx>y 

i_ f h R{y) dxdy~- [ A dx) \a + ) + 0(t 4 

T J x~>y ^ J 
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where the integrations are from e to A, subject to the path ordering constraints as indicated. 



Order t 

Our strategy will be to expand R(x) in modes around zero and analyse which modes give divergent 
contributions to the integral. To do so, first express R(x) in the AT-basis (I2.19[) . 

R(x) = K 3 (x) \(k+{x) K-(xj) A . (A.10) 

The field K 3 (x) has integer modes and K ± (x) half-integer modes. We decompose 

K*{x) = K%{x)+K%{x), Kl(x)= ^ m ~ lK3 rm K K X ) = E *~ m " 1 *m, 

m£Z> mGZ <0 

K ± (x) = K±{x)+K±(x), K±(x)= x~ r ~ l Kf, K±{x) = ^ x~ r - l Kf . 

(A.ll) 

We also set R<{x) and R>(x) to be be given by IjA.lOf) with all fields replaced by their <-part, 
respectively their >- or >-part; the f- term is part of R>{x). The coefficient of t in l|A.9|l becomes 

~ J R(x)dx\a+) =-jR < (x)dx\a + ) + -J (K 3 - \fx)x- l dx \a+) . (A.12) 

The integral over R < (x) involves only powers x r with r > — | and has a finite e — > limit. The 
second integral has a log-divergence for e — > which has to be cancelled. This is achieved by setting 
fx = 1. With this choice for fx we have 

R>(x)\<r+) =Q. (A.13) 

Order t 2 

For the second order computation, we need to know the commutators [K>(x), K b {y)] for x > y. 
These can be computed from the commutation relations of the A^-modes, which are just the same 
as those of the J^-modes given in Q2.14[l . One finds, for v = ±, 



[KUx),K\y)\ 



1 



2{x - yf 
v 



[Kl{x),K»{y)} = K»(y) 

x — y 

y v 



[K»(x),K 3 (y)] = - J y - K"{y) 

V x x — y 



(A.14) 



For the coefficient of t 2 in 1A.9I) wc then obtain 



dy - - [ dx ) |er_, 



(4 / R(x)R(y)dxi, a 

= (-^/ R < {x)R < {y)dxdy+\ ( [R>(x), R(y)] dxdy - A(ln A - lne)) k + ) 
J x>y ir 2 J x>y L 4tt / 

i / Jr*W + A-M i / i (A . 15) 



+ ^ In e+ (finite e^0))|er + ) 
= f J_lne+ ^ln e+ (finite e-fO)W; 

V 47T^ 47T / 
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where in step 1) we replaced R(x) by R<(x) + R>(x) and used (|A.13(1 . For step 2) note that 
the integral over R < (x)R < (y) has a finite e — > limit. The commutator [R>(x), R(y)] has been 
evaluated with the help of 1A.14J) . For step 3), the singular contributions to the resulting integrals 
have to be extracted. The most singular term in the first integrand (expand K u (y) in modes) is 
(y/xy(y/x + ^/y))^ 1 = 4d x dy((y/x + <Jy) ln(y/x + ^/y)), so that the integral is regular for e^O. The 
integrand of the second integral is —Ad x d y ln(i/a; + ^/y) so that there is a In e singularity. In order to 
have a finite s^O limit at order t 2 we need to choose f% = —n^ 1 . 



Order t 3 

The coefficient of t 3 in (|A~3jl can be written in the form A\ + A 2 + A3 with 

Al = ~3 [ R{x)R{y)R{z) dx dy dz \a+) 

1 x>y>z 



TV 

A - h 

47T 2 



rm + m\ dxdyW+) (A.i6) 

V v x J 



x>y V V 



The singular contributions to A 2 and A 3 are easy to evaluate, 

A 2 = ]n(e)-Pn I R<(x) dx \a+) + (finite e^O) , A 3 = ln(e)^-\a + ) + (finite e^O) . (A.17) 
Air z J Air 

Extracting the singular part of A\ is some work. One first uses the commutators l|A.14(l to remove all 
positive mode parts of the fields K a . One obtains a sum of terms, with each term being a function 
in x, y, z multiplying a state of the form, for u, v € x, y, z, 

\a+), K%(u)\a+), K a < (u)K b < (v)\a + ) , K a < (u)K b < (v)K c < (u)\a + ) . (A.18) 

It turns out that only the coefficients of the first two terms lead to singular behaviour as e — > 0. In 
fact on finds 

M = -htt f J0n$w dx dy dz + (finitc £ ^ 0) ' (A - 19) 

where 

if dx dy dz ln( ) + (fi 't 0) (A 20) 

J x>y>z 

Altogether 

Ai + A 2 + A 3 

= ln(e) (~^ + 4hl R<{x)dx+ i^l R<ix)dx+ %%) |fT+> + (finite e ^ 0) ' ! A " ' 

which has a finite limit iff f 2 = — 7r _1 and / 3 = i. Note that the required value for f 2 agrees with 
the one obtained in the order t 2 computation. 
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